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Abstract—This article reviews some of the research 
carried out by our group on modeling of metamaterials 
whose building blocks are not very tightly coupled. To 
solve the problems involving complex metamaterials, 
which may either be periodic with slowly convergent 
behaviors or large non-periodic arrays, we apply the 
Characteristics Basis Function Method (CBFM) which is 
a powerful paradigm for handling such large complex 
problems. The problems are formulated by using integral 
equations and the Method of Moments (MoM) in 
conjunction with the Characteristic Basis Function 
Method (CBFM) to significantly reduce the number of 
unknowns. We apply the technique to solve problems 
involving different plasmonic metamaterials including 
periodic nanorods, finite random arrays of rods as well as 
a structure which has a general loop geometry building 
block. 


Index Terms— Metamaterials, 
characteristic basis functions 


plasmonics, 


I. INTRODUCTION 

Arrays of plasmonic antennas exhibit 
interesting physical attributes and find a wide 
range of applications in photonics. Light- 
trapping arrays of plasmonic nanoparticles 
that can increase the photo-current spectral 
response of thin film solar cells have been 
studied systematically in [l1]. Arrays of 
hetero-structure gold nanorods loaded with 
poly  (3-hexylthiophene) (P3HT), an 
electronically-active semiconductor [2], can 
concentrate the light in the nanoscale gap 
between the two antenna arms [3]. One- 
dimensional arrays of plasmonic nanoparticles 
have been reported to operate as optical 
waveguides supporting different guided 


modes with different polarization properties 
[4]. An array of rectangular gold nanorod 
antennas surrounded by an insulating S10, 
region and deposited on a semiconductor 
interface can act like a Schottky diode to 
generate hot electron-hole pairs, resulting in a 
photocurrent [5]. It has been demonstrated in 
[6] that nano-optical Yagi-Uda antennas, 
comprising of plasmonic nanorods, can both 
enhance and direct the interaction of a single 
quantum emitter with electromagnetic fields. 
Also, arrays of dielectric-plasmonic core-shell 
nanoparticles can be used to shape the near- 
field pattern, enabling sub-wavelength optical 
focusing [7]. Another interesting example 
worth mentioning is the immobilization of 
bacteria at designated positions in space to 
study their metabolism without damaging 
them by utilizing an array of plasmonic 
nanoantennas that produces strong gradients 
of light intensity for the trapping mechanism 


[8]. 


While these structures are interesting 
because of their physics, simulating them is 
very challenging because of their large size 
since they are not necessarily periodic. The art 
of metamaterial design is based on judicious 
choices of periodic or truncated periodic type 
with elements that are sub-wavelength-sized 
structures consisting of various materials and 
geometries often regarded as unit cells or 
building blocks. Obviously the design relies 
on our understanding of the physics involved 
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in the function of the building blocks. The 
design process usually begins with the 
analysis of isolated blocks and proceeds with 
the modified configurations to achieve novel 
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Fig. 1: (a) Two-Dimensional array of plasmonic nanorods illuminated by an oblique incident plane wave. The 
periods along x- and y-axis are A, and A,, (b) The plasmonic nanorod made of a Drude material as the unit-cell, (c) 
Real part and imaginary part of the permittivity of silver over a desired optical frequency range. 


physics. While realizing equivalent bulk 
properties, such as doubly-negative material 
indices, 1s often our main goal, such properties 
can only be guaranteed through careful design 
and analysis of the individual building blocks 
as well their periodic or quasi-periodic 
versions. More recently, large arrays of 
metamaterials are finding interesting 
applications on platforms such as solar cells, 
aircraft antennas, antenna radomes and optical 
sensors. This further enhances the need for 
efficient design tools for handling such large 
truncated metamaterial arrays, and the 
availability of suitable computational tools is 
of paramount importance to the successful 
design and analysis of metamaterials. While a 
large variety of computational methods exist 
that can potentially address the requirements 
of metamaterial design, in this work we focus 
on the Method of Moments (MoM) and a 
projection-based acceleration of MoM. 


In the characteristics basis function 
paradigm, an element is first isolated and 
solved for plane wave excitations incident 
from different angles. Next, a singular value 
decomposition (SVD) is applied to extract a 
few basis functions that capture the physics of 
the element, reducing the number of 


unknowns significantly (sometimes by orders 
of magnitude) in the process. Next, the 
solution to the array problem is expressed in 
terms of the CBFs for the elements, and their 
coefficients are found. 


In this work we employ the CBFM to model 
three different configurations: (1) a periodic 
array of plasmonic antennas; (11) random and 
Fibonacci array of plasmonic rods; and (iii) a 
metamaterial comprising of plasmonic rings. 


Il. PERIODIC ARRAY OF 
PLASMONIC RODS 


Figure 1 illustrates the geometry of a 
periodic array of silver nanorods with a=7.5 
nm, L=150 nm, and periodicities A,=150 nm 
and A,=250 nm. The nanorods are orientated 
along the y-axis and are illuminated by a plane 
wave with parallel polarization, as shown in 
Fig. 2(a). As depicted in Fig. 2(b), it can be 
shown for an electrically small dielectric 
sphere, with an arbitrary permittivity €, the 
scattering characteristics can be derived from 
an equivalent electric dipole moment 


I = 4aj(k,a) (e, -1)E, /[nk; (£, +2)] , where 
ko and y are the free-space wave number and 
intrinsic impedance, respectively, and KE, is 


the incident electric field measured at the 
center of the sphere. We utilize the same 
concept for modeling longitudinal polarization 
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authors have shown that by discretizing the 
current using only a few piecewise sinusoidal 


Electric Dipole 


Du- (ka 

















é.—l 
E 
ie ° 


E, 
Z Polarization Current 
k inc H 





(d) 


Fig. 2: (a) Characterization of a single nanorod illuminated by a Parallel polarization, (b) Modeling the scattering 
problem of plasmonic nanorod originated from Dipole Mode (DM) approach, (c) Discritization of the polarization 
current by piecewise sinusoidal macro basis functions, (d) Full triangular macro basis function including right and 


left half-MBFs. 


Macro Basis Functions (MBFs), as depicted 
in Fig. 2(c), and by employing Galerkin 
testing on the surface of the nanorod, one can 
derive an accurate, fast and singularity free 
computational paradigm to solve the problem 
at hand and to obtain the unknown current. 
The current induced on the axis of nanorod 
can be derived by solving the matrix equation 
ere ye av O 


res 


where Teas lis a column matrix that includes 
the coefficients of each MBFs, and 


Ym+ 


[Ze In = ftn O), dy, (2) 


Ym-1 


Ym+1 


VE n= [E in ay, (3) 


Ym-1 


and = t(yv)=sin(ko(A-|y-vn|)) is the mth 
piecewise sinusoidal MBF. As depicted in Fig. 
2(c), the first and last MBFs are one-half while 
the rest are full piecewise sinusoidal functions. 
The length of the domain of the half MBF is 
H, and y,=-L/2+(m-1)H is the location of cusp 
of each MBF. The electric field radiated by 


MBF#n is denoted by G” in (2) and is 
defined as: 


GPE Right (y y ) n=1 1) 
Ge = G” AI (y —y,) n = N wer > 
GMBF Ful (y _ y,) Otherwise (4) 


where G's" is the electric field, radiated 


by the right half triangular-sinusoidal depicted 
in Fig. 2(d). Solution of (1) enables us to 
obtain the unknown current on the nanorod, 
and we can solve for this current for various 
incident angles and frequencies (see Fig. 3). 


Our next task is to account for the mutual 
interactions between the elements and obtain 
the polarization current for the elements 
within the array. Toward this end, we use a set 
of high-level entire-domain basis functions 
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Fig. 3: Maximum of the magnitude of the longitudinal 


current versus the frequency and the azimuthal incident 
angle (o) for the excitation shown in Fig. 2(a). 


called the Characteristic Basis Functions 
(CBFs) . To obtain the CBFs for a scatterer, 
we illuminate an object as a set of plane waves 
[PWi, PWo,...,PW,, ]. Let us assume that J,( 


r') is x-, y- or z-component of the induced 
current vector within the rod, for a specific 
plane wave, e.g., #n, PW,. Next, we construct 
the column matrix [Jn] which contains N; 
samples of J,(1r') at N, different points within 
the scatterer, r, , Isis N. Then we construct 


a rectangle matrix [J] whose n column is 
[Ja]: 


[J] =[F Lab. y, I (5) 


The CBFs can be obtained by using the 
Singular Value Decomposition (SVD) 
factorization [10] of [J]. Figure 4 illustrates 
the magnitudes, phases and corresponding 
singular values of the CBFs for the plasmonic 
nanorod at f=260 THz (first resonance), /=400 
THz, and /=440 THz (second resonance). The 
criterion for retaining the CBFs is based on the 
threshold for the normalized singular value, 
which we set to be above -30 dB, Le., 
10xlog(Snn/S11)>-30. It can be observed that at 
f=260 THz, only two CBFs are needed, 


whereas at higher frequencies, /=400 and 440 
THz, three CBFs are needed. The first CBF, 
CBF#1, at /=260 and 400 THz, has one peak, 
while it possesses two peaks at f—-440 THz. 
We also observe that the higher-order CBFs 
have more variations in both magnitude and 
phase compared to the lower order CBFs. 

To solve for the periodic array with 
unknown current /(y), we need to impose the 
polarization current condition [9] on the center 
element: 


jae,(€, -l(E" +E") =I(y)/ (ra) (©) 
-L/2sysL/2 


where Ee is the y-component of the scattered 


electric field, generated by the infinite 
elements of the array. To solve (6), we express 


I(y)in terms of the CBFs: 
Nc 
CBF 
y= XL hO), (7) 
n=] 


where Nc is the total number of CBFs. This 
number is relatively small, for instance Nc=3 
at f-400 THz when we impose the -30 dB 
criterion on the associated singular values. The 
n™ CBF is denoted by f(y) in (7), and 1" is 
its unknown coefficient. The scattered electric 
field, Ee , can be expressed as an infinite 
summation of the fields radiated by each of 
the Nc CBFs residing on each of the elements, 
including the center element. By taking into 
account the inter-element phasing factor ep, 
introduced in (5), we can write EY“ as: 


Nc +00 
scat CBF CBF ,n 
Ee = Ee eGo (8) 
n=l P.q=-© 


where Go ” is the y-component of electric 


field radiated by the nth CBF on the (p,qg)th 
element. 


Since each of the CBFs is a superposition of 
MBEFs with known coefficients, we can derive 































linear combination of expressions for the 
fields radiated by the full- and half-MBFs that 
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Fig. 4: Magnitude and phase of CBFs and singular values for plasmonic nanorod at (a) f=260 THz, (b) f=400 THz, and 
(c) f=440 THz. 


replacing (7) and (8) in (6), we test the 
equation (6), by using m'th CBF, lsmsJN_, 
at the center element (x,y,z)=(0,y,a). Such a 


Galerkin’s testing [11] leads to an NcxNc 
matrix equation which reads: 


(aa B [cP | |r |- [vo] (9) 
where 
CBF B f < CBF ,n 
Ze bin = f Sn(Y) denon Ws — (10) 
-L/2 P,q=-* 
(a ea fS DLOdy, (11) 
res dmn ma k, (E, ~1) m n 


-L/2 


Li2 


[VO], = SE OSO, 


-L/2 


(12) 


and | is the column vector representing 
the Nc unknown coefficients [°* . The 


algorithm we employ to evaluate the infinite 
series in (10), comprises of some successt?? 
stages using the ring concept described above. 
At stage #M we need to consider the truncated 
(2M-1)x(2M-1) array including the rings #1, 
#2,.... #M to find the corresponded coefficients 


rc po. Thus, at stage #M, we solve the 
matrix equation (9) by using (10-12), except 
that instead of using a doubly-infinite series to 


CBF 


construct [zę |in (10), we use a truncated 


series in the spatial domain as follows: 


L/2 
PRD f L(y) Se Gay, (13) 
L2 p.q=-M 


Since the elements of rc r go up and 
down synchronously, we increase M, the 
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Fig. 5: Reflection and transmission coefficients (in 
percentage) for normal and oblique incidence cases 
compared with FDTD. The FDTD result is with limited 
accuracy due to the mesh-size. (a) Normal incidence, 
(b) Oblique incidence. 


number of rings, as indicated above, the shape 
of the composite CBF does not change as the 
number of rings is increased. Thus, we only 
need to determine its weight coefficient. 

The reflection and transmission coefficients 
for both normal and oblique excitation cases 
are plotted in Fig. 5, and are compared with 
full-wave FDTD simulations. For the FDTD 
modeling, a fine mesh must be used to model 
the thin cross-section to achieve accurate 
results. Considering the large aspect ratio of 


the rod and its dispersive material property, 
the simulation is time consuming and the 
FDTD simulation takes approximately 10 
hours. A coarser mesh provides a faster 
computation (of about 5 hours) but is less 
accurate. As observed in Fig. 5(a), our CBF 
model provides a good comparison with the 
FDTD fine mesh results. For oblique 
excitation, we use Sin/Cos single frequency 
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Cc) 
Fig. 6: Reflection and transmission coefficients (in 
percentage) for different incident angles and different 
polarizations, (a) For three different values of 0 =25", 
55°, and 85°. The values of ọ=0 and w=0 are fixed, (b) 
For three different values of ¢=30°, 90°, and 120°. The 
values of 0=60° and w=0" are fixed, (c) For three values 


of w=30°, 100°, and 120°. The values of 0=30° and 
ọ=60° are fixed. 


calculation [12] in FDTD with coarse 
meshing, and the result is shown in Fig. 5(b). 
This calculation, performed at 24 points, 
requires approximately 5 days to complete. In 
comparison, the CBFM is considerably faster 
by orders of magnitude, and it takes only 
about 30 seconds for 73 samples of this 


























ME eh o 
1000H--H--HH-HH-H-H INES 
HEEL EEE ETE) £o IZANTA © 
HET ET ETE LIVAALIZUANNZ 
500H--H--HH-HH-H-H g 5OOHTLHINAMAVLH! k 
= HETET ETE = ATTA 
€ PEEVE = HIATININTL 
TTT) E aa 
> HPETET ETT > VATLULAN LY 
-500H HEHEHE -500dAH AAA 
HL ALA 
1000 EEE EEE EEE EA AOOOH EET AIA 
PTTT ETE UI TTT FALLAN 
-500 0 500 -500 0 500 
x [nm] x [nm] 
Ca) (b) 


FIG. 7: Finite array of plasmonic nanorods (a) periodic 
array, (b) disordered array. 


frequency to generate an accurate result. It is 
evident that compared to the FDTD, the 
CBFM is relatively independent of the angle 
of incidence, due mainly to the fact that the 
CBFs include all the features of different 
excitations. With regard to the convergence of 
spectral series in (13), the truncation number 
equal to N7=20 or 25 is sufficient regardless of 
the angle of incidence, because of the 
decaying exponential for higher order Floquet 
modes and the decaying behavior of the 
Fourier Transform of the CBFs for relatively 
large |ķ%|. For both normal and oblique 
incidence cases, the resonance frequency 
occurs at f= 262 THz, which is close to the 
resonance of the single nanorod f=260 THz. 
This is mainly because the couplings between 
the elements are not strong. At resonance, the 
reflection coefficients for the normal and 
oblique incidence cases are 76% and 54%, 
respectively, whereas the corresponding 
transmission coefficients are 2% and 29%. 
Hence, the absorption loss in the normal case 
is 22%, and 17% in the oblique case. 


Additional results for other oblique incidence 
cases are depicted in Fig. 6 for different 
combinations of y, 9, and ọ. 


HI. RANDOM ARRAY OF 
NANORODS 


In this section we address the problem of 
modeling a finite random array of plasmonic 
elements. Fig. 7(a) demonstrates the geometry 
of a 15x15 array of plasmonic nanorods. The 
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FIG. 8: Optical scattering performance for arrays of 
plasmonic nanorods for finite periodic, disordered and 
doubly infinite configurations. In case I, A,=100 nm 


and A,=170nm, and in case 2, A,= Ay =250 nm. 


unit cell, containing a plasmonic nanorod, 1s 
rectangular in shape with dimensions of 
A,=100 nm and A,=170 nm. The plasmonic 
nanorod has the dimensions of Z=150 nm 
length and radius a=7.5 nm. The rod is made 
of silver, whose permittivity can be 
characterized with the Drude model, 1.e.Ć, 
£, =E- fy US - itd]. with en =5. f 
=2175 THz and fa =4.35 THz. The excitation 
is a plane wave illuminating the nanorod at 
broadside and its polarization is vertical (along 
the y-axis). Fig. 7(b) shows the same array 
with a random disorder in which each of the 
nanorods is randomly tilted around its 
centroid. The tilt angle between the directions 
of each of the nanorods and the y-axis varies 
randomly in the range of +30°. Such a disorder 
in the geometry can be attributed to 
imperfections in the fabrication or the non- 
periodicity and may be introduced deliberately 
to investigate non-periodic phenomenon. 


Once again, we apply the CBFM proposed 
in the previous section to solve such a large 
array problem. Fig. 8 illustrates the magnitude 
of the scattered electric field, |E], observed at 
r = 3| um]z (2.5 at 250THz) for the periodic 


and disordered finite arrays (15x15) for two 
different cases. In case I, the inter-element 
spacing is the same as mentioned before, but 
in case II we increase the spacing to A= A, 
=250 nm. The first thing we observe is that the 
coupling increases as we widen the bandwidth. 
In general, Case I has a wider bandwidth than 
does Case II. It 1s also interesting to note that 
the resonant frequency for the array with 
closer elements is lower than that of the array 
with increased element spacing. The resonant 
frequency of the periodic array in case II is 
close to the frequency at which the length of 
the nanorod is one-half of the effective 
wavelength, Ap /2 [13]. The +30° disorder in 


the alignments of the nanorods affects the 
performance of the Case II only slightly since 
the coupling levels are less in comparison to 
those in Case I (less interaction). For Case I, 
the disorder can shift the resonance frequency 
upwardly for which the |E| is maximum. For 
the sake of comparison we also present the 
results for the case of an infinite periodic array 
obtained by using the closed form formula for 
the reflection coefficient in. 

Next we solve the problem of a 2D 
Fibonacci array [14, 15] by using the CBFM. 
In one spatial direction, 1.e., x- or y-direction, a 
Fibonacci array of plasmonic nanorod 
antennas is generated by using binary 
Fibonacci sequence consisting of two seed 
letters A and B, where A represents one 
plasmonic nanorod with the length L4 and B 
represents the other plasmonic nanorod with a 
different length, Lz. The Fibonacci sequence 
of order n>2, Fn, can be generated by using the 
recursive formula, F,=F,-1 Fy-2, where F7=A, 
and Fo=B. For example, the one-dimensional 
Fibonacci sequence of order 6 has the form: 


F, = ABAABABAABAAB. (14) 


The length of the sequence Fẹ is 14. 


We design two plasmonic nanorods A and 
B, such that the resonant frequencies are at f4 
=260 THz and fg=1.05* f4 =273 THz. Such a 
design can be carried out approximately by 
taking advantage of the algebraic closed-form 
formula for the effective wavelength of a 
plasmonic nanorod in [13] 
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FIG. 9: The magnitude of the scattered field by the 


designed A and B plasmonic nanorods to get the 
resonance at f4=260 THz and f/g=1.05 f4=273 THz. 
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FIG. 10: Fibonacci quasi-lattice of order 6, 13x13. The 
array comprises of two different plasmonic nanorods, 
blue nanorod with lower resonant and red with higher 
resonant frequencies. 


where the physical length of the nanorod is 
one-half of the effective wavelength. 
However, to increase the accuracy of the 
design, we use the CBFM to find the 
frequency response of the magnitude of the 
scattered electric field by the plasmonic 


nanorod illuminated at broadside by a 
vertically polarized plane wave. The 
maximum of the scattered field over an optical 
range of frequency occurs at the resonant 
frequency. For the same radius a=7.5 nm and 
the same material, silver, for the both A and B 
nanorods we obtain L4 =150 nm and Lg =139 
nm, for the previously mentioned desired 
resonant frequencies (see Fig. 9). The 
observation point is at r =3[um]z, assuming 


that the nanorod is located along the y-axis. 
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FIG. 11: The optical performance of the Fibonacci and 
alternative quasi-lattices of plasmonic nanorods. 
Bringing the elements closer to each other can merge 
the two resonances. 


Having designed the two nanorods A and 
B, we can next construct the 13x13 Fibonacci 
quasi-lattice of order 6, demonstrated in Fig. 
10. The arrangement of the A and B nanorods 
located in the first row and column of the 
array in Fig. 10 is according to Fs in (2). For 
other columns and rows, the arrangement is 
either Fs or F, (the complement of Fs for 


which A and B are simply interchanged), 
depending on the first element (A or B) with 
which the column or row starts. 

The optical performance of the 13x13 
Fibonacci quasi-lattice of nanorods is 
illustrated in Fig. 11 by using the CBFM, 
described before. The result is also compared 
with the 13x13 alternative array which instead 
of Fs it uses the alternative sequence S 
=ABABABABABABA, of the length 13, and 


its complement S . The inter-element spacing 


in both the x-and y-directions is Ax= Ay =200 
nm. For both the Fibonacci and the alternative 
arrays, the frequency response can be 
described as a blend of two resonant responses 
illustrated in Fig. 9. For the Fibonacci case, 
the second resonance has a slightly wider band 
and it is closer to the first resonance. 
Moreover, in Fig. 11, we characterize the 
optical performance of the 53x13 alternative 
array, in which we bring the elements closer 
along the x-direction by choosing A,=50 and 
Ay =200 nm, but we increase the number of 
elements in x-direction to have approximately 
the same aperture size as the 13x13 alternative 
array. We see that the two resonances are 
merged and a wider bandwidth is achieved 
when we bring the elements closer to each 
other. 

We need to stress, once again, that our 
model can characterize the performance of the 
finite arrays (disorder or Fibonacci) extremely 
quickly and accurately. The modeling of the 
large arrays presented herein takes about 6 
min/40 frequency samples on the machine 
Alienware 17-3.47 GHz CPU with 8.00 GB 
RAM over the entire frequency band. Our 
technique can be used as a powerful engine for 
efficient analysis and novel design 
optimization of large and  non-periodic 
plasmonic configurations. 


IV. MORE COMPLEX 
CONFIGURATIONS 


For a more complex geometry, which is a 
building block for a metamaterial structure, we 
again employ an integral equation method 
based on MoM, to analyze the building block 
and extract the CBF modes. A variety of 
integral equation (IE) formulations exist in the 
literature for penetrable objects. Most 
formulations rely on the well-known Stratton- 
Chu type of formulas for the field 
representation formulas [16]. The results 
reported in this article have been obtained 
using the JMCFIE formulation as it is 
reflected in JMCFIE [17, 18]. 
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Fig. 12: Plasmonic ring. The frequency response of the 
surface field energy on the 40x40 array. The È. and 
K, in the legend 


propagation direction of the incident field. The arrow 
pointing to the peak of the curve indicates the actual 


signify the polarization and 


value of the plasmonic U around the resonance 


frequency. 


In JMCFIE, J and M represent the 
equivalent electrical and magnetic surface 
currents on the material interface separating 


two homogeneous material regions while £” 


and H” represent the incident fields in the 


region 7. At the same time, L and K are the 
integral operators used in the field 
representation formulas and are defined in the 
EFIO (electric field integral operator) and 
MFIO (magnetic field integral operator), 
respectively, as follows. 
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Fig. 13: Plasmonic ring. The real part of the tangential 
field due to the first six MBFs generated for the array at 


f =S50THz. The E, and the H, captions signify 


electrical and magnetic field representation of the i” 
MBF 


K()={ Vig(ir-r' pJ(r)dS’ (17) 


We apply the model to solve the problem of 
a metamaterial structure comprising of 
plasmonic rings. Its complex configuration 
presents a greater challenge from the point of 
view of computational modeling but provides 
unique features to control resonances. It has 
been shown that the behavior of the loop 
structure can be successfully controlled by 
controlling the aspect ratios associated with 
the geometry in the transverse and vertical 
planes [19]. The plasmonic ring discussed here 
is a ring with internal radius of 20nm , 
external radius of 40nm and thickness of 
30nm . The lattice constant of the array 
structure is chosen to be equal to 300nm. As 
depicted in Fig. 12, with the excitations given 
as E” =p "Y and E" =e *® , the 
structure shows resonances at around 540THz 
and SSSTHz, respectively. 


The first six macrobasis functions are listed 
in Fig. 13. The first two modes have identical 
patterns but differ from each other by a 90 
rotation around the Z axis. The second two 
modes are also identical in pattern but differ 
by a 90 rotation around the 2 axis. 
Nevertheless, the third and the fourth mode, 
respectively, represent a magnetic and an 
electric dipole mode along 2. 


Conclusion 


In this work, we have reviewed the concept of 
Characteristics Basis Function Method 
(CBFM) for solving the problem of 
metamaterials that are large finite arrays. The 
fundamental concept upon which CBFM is 
based, 1s to extract the mode profile for the 
building blocks. We isolate a building block 
and excite with plane waves from different 
angles, and down-select a few modes by using 


the SVD to capture the physics of that block. 
This method is powerful as well as efficient 
both in terms of computational time and 
memory. The CBFs need to be generated first, 
and the reduced impedance matrix is formed 
next. Both of the procedures may be time- 
consuming when compared to the 
conventional MoM if the array is not 
sufficiently large, though dividing the problem 
into smaller sub-problems can help with this 
issue. The number of blocks affect the 
generation of CBFs as well as the filling of the 
reduced impedance matrix. Understanding 
some of the issues 1s helpful for improving the 
performance of the CBFM algorithm even 
further. 
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